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SUMMARY 


In  this  paper  ie  is  shown  Chat  some  distributions  of  the  matrix 
variates  and  latent  roots  arising  in  the  multivariate  Behrens- 
Fisher  discriminant  problem  can  be  explicitly  expressed  in  terms 
of  the  invariant  polynomials  with  two  matrix  arguments,  due  to 
A.  W.  Davis,  extending  the  zonal  polynomials  of  matrix  argument. 
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1 .  INTRODUCTION 
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We  are  interested  in  the  discrimination  problem  of  two  /  '  '  * 

multivariate  normal  populations  under  the  heterogeneity  of  pop-  i  . 

ulation  covariance  matrices.  Distributional  Problems  for  the 
univariate  case  have  been  investigated  by  many  authors  (e.g., 

McCullough,  Gurland  and  Rosenberg  (I960),  Ray  and  Pitman  (1961) 
and  Welch  (1947) ) ,  and  Gurland  and  McCullough  (1962)  and  Wehrhahn 
and  Ogawa  (1975)  took  into  consideration  the  preliminary  testing 
procedures  for  the  null  hypothesis  of  the  equality  of  the  variances. 

Yao  (1965)  considered  approximate  distributions  of  a  Hotelling's 

2 

T  -type  statistic  in  the  multivariate  case.  However,  to  the 
author's  knowledge  very  little  work  has  been  published  on  the  multivariate 
Behrens-Fisher  discriminant  problem.  Our  Behrens-Fisher  matrix  statis¬ 
tic  for  the  multivariate  case  may  be  in  the  form 


D  -  (Sx  +  S2)"lsS0(S1  +  S2)”** 


where  the  m  *  m  matrices  are  independently  distributed  as,  in 


(1.1) 
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general,  noncentral  Wlshart  Wm(n^,  1^,  D^) ,  1  «  0,  1,  2.  In 
practical  terms,  Sj  and  S2  are  sample  error  ss.  matrices,  pro¬ 
portional  to  estimates  of  the  distinct  population  covariance 
matrices,  and  Sq  Is  an  effects  ss.  matrix,  estimating  tho 

difference  of  the  population  mean  matrices.  In  the  case  when 
2  2 

Ei  "  ^2  “  o  E0 (o  >  0)  and  Qi  *  H2  “  0,  D  is  the  well-known  MANOVA 

matrix  (central  if  flo  *  0  and  noncentral  if  ftp  0) .  Hence,  D 

is  an  'extension'  of  the  MANOVA  matrix  in  discriminant  analysis. 

-Vt  -W 

The  matrix  A  «*  (Sj  +  S2)  Sj (Sj  +  S2)  is  a  suitable  preliminary 
test  matrix  for  the  hypothesis  «  E2,  where  in  this  paper  we  assume 
n,  -  02  -  0  whenever  the  preliminary  testing  procedures  are  considered. 
When  Ej  $  E2,  D  and  A  are  dependent,  and  we  may  be  interested  in  the 
conditional  or  unconditional  distributions  of  the  roots,  of  D,  or 
suitable  functions  of  them  (conditional  on  the  roots  of  A). 

In  this  paper,  it  is  shown  that  some  distributions  of  the  matrix 
variates  and  latent  roots  arising  in  the  multivariate  Behrens-Fisher 
discriminant  problem,  for  the  null  case.  Sip  «  0,  can  be  explicitly 
expressed  in  terms  of  invariant  polynomials  with  two  matrix  arguments. 
The  invariant  polynomials  with  two  matrix  arguments  have  been  defined  by 
Davis  (1979a),  (1979b),  in  the  need  of  extending  the  zonal  polynomials 
due  to  Constantine  (1963)  and  James  (1964).  Properties,  relationships 
between  them  and  applications  in  multivariate  distribution  theory  are 
discussed  in  the  same  papers  of  Davis.  It  is  noted  here  that  explicit 


forms  for  those  distributions  for  the  non-null  case,  Hq  ^  l»e*» 

when  population  mean  matrices  are  different,  can  be  obtained  by  in¬ 
troducing  invariant  polynomials  with  three  matrix  arguments,  extending 
the  work  of  Davis.  The  investigation  of  invariant  polynomials  with 
larger  number  of  matrix  arguments  will  be  presented  in  subsequent  papers 
(see  e.g.,  Chikuse  (1980)). 

The  invariant  polynomials  with  two  matrix  arguments  are  briefly 
discussed  in  Section  2.  The  joint  distributions  of  the  roots  of 
D  and  A  are  derived  in  Section  3,  and  the  'marginal’  (and  conditional) 
distributions  of  the  roots  of  D  (conditional  on  those  of  A)  are 
investigated  in  Section  4,  for  the  null  case,  ftp  "  °* 
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2.  INVARIANT  POLYNOMIALS  WITH  TWO  MATRIX  ARGUMENTS 

Davis  (1979a),  (1979b)  has  defined  a  class  of  homogeneous 
ic  X 

polynomials  C  *  (X,  Y)  of  degrees  k  and  l  in  the  elements  of  the 

<0 

m  x  m  symmetric  complex  matrices  X  and  Y,  invariant  under  the  simultaneous 
transformations 

X  -*•  H'XH,  Y  •+  H'YH,  H  e  0(m),  (2.1) 

where  0(m)  is  the  group  of  m  *  m  orthogonal  matrices.  These  satisfy 
the  basic  relationship 

(  C  (AH'Xri)C. (BH'YH)dH  -  L  C^,X(A,  B)C^,A(X,  Y)/C.(I),  (2.2) 

J0(m)  K  A  ♦eic*  *  'P  ♦  + 

where  C^,  C^,  are  zonal  polynomials,  indexed  by  the  ordered  partitions 

k,  A,  if>  of  the  nonnegative  integers  k,  i,  f  «*  k  +  SL  respectively  Into  not 

more  than  m  parts  (for  the  zonal  polynomials,  see  e.g.,  Constantice  (1963) 

and  James  (1964)).  Letting  G£(m,  R)  denote  the  group  of  m  *  m  real  non- 

sipgular  matrices,  * e  <«A'  signifies  that  the  irreducible  representation 

of  G!(m,  R)  indexed  by  2<j>  occurs  in  the  decomposition  of  the  Kronecker 

product  2k  0  2A  of  the  irreducible  representations  indexed  by  2k,  2X.  The 

K  X 

properties  and  relationships  satisfied  by  the  C  *  ,  which  are  especially 
utilized  in  the  later  sections,  are  summarized  in  the  following  (see 
Davis  (1979a),  (1979b)  for  proofs  and  details): 

CJ*X(X,  X)  =  0*’XC^),  where  0*’A«  C*’A(I,  I)/C^(I),  (2.3) 

C*,A(X,  I)  *>  [9^’AC^(I)/Ck(I)]Ck(X),  and  similarly  for  C*’A(I,  X),  (2.4) 

C*»°(X,  Y)d*f  Ck(X),  CX°,A(X,  Y)d=f  CX(Y),  (2.5) 
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ir  . 


C*’X(aX,  @Y) 
4> 

CK(X)CX(X)  - 
where  g*>x  -  I*. 


5  a^S^C^’^CX,  Y)  (a,  3  complex  constants), 

vx)  ■  W^v*- 

r  i  ?- 

(8**  )  »  which  is  defined  in  Constantine  (1966, 


Eq.  (27)), 

f  C*‘X(AH’XH,  AH' YH) dll  -  C*,X(X,  Y)C  (A)/C  (I)  , 

J0(m)  ♦  ♦  *  * 

I  etr  (^R)|R|a"PcJ‘X(XR,  YR)dR  -  ^(a,  ♦)  |W|  "*cJ*X(Xtf"1,  YW*1), 


(2.6) 

(2.7) 


(2.8) 

(2.  9 


where  p  ■  (m  +  l)/2, 

{  C*?*X(A'HX1IA,  B)dH  =  C*f  ’ X(A* A,  B)C  (X)/C  (I) ,  (2.10) 

JOC m)  *  * 

|glR|t~P|I  -  R|U  PciCAR)<iR  -  rm(p)|x|  l~„  krm(t.« 

(-u  +  p)  G*’V’X(X,  AX)/k!r  (t+p,  <!>),  (2.11) 

<  (p  <p  in 

[I|R|t_P|i  -  RfV-V  i-R)dR  -  rm(t,  »)r,(«,  x>ej«\«>/rm(t+u,  «>. 

0  (2.12) 

A  multivariate  generalization  L^(X,  A)  of  the  Laguerre  polynomial  due 

to  Khatri  (19.77)  has  an  expansion 

>><*.  *>  •  +  >>,  L,X(*«.X)(k>ed,Xc?X(A'  •itA)/(t  +  RV  (2‘13) 

The  following  lerrana  is  also  useful. 

LEMMA  2.1. 


<1> 


»<R+T<I 


a  -p  b-p  c  -p  . 

I  -  R  -  T|  |R|  |T|  cj»  (R,  T)dRdT 

-  r  (a)T  (I,  <)r  (c,  A)e*’xc.(i)/r  (a*t>+c,  *), 

m  m  m  <P9  ra  ’  T,» 


(2.14) 
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rX  a-f>  k,X 

(ii)  Jo  |r|  (ara*  ,  a)  dR  =  rni(P)rm(a,K)|x|ac^'A(AXA',B)/rn(a+p,  k). 


(2.15) 


Proof.  (1)  We  evaluate  the  Integral  form 


r-f  [  f  etr-(X  +  Y  +  Z)|x|a"P|Yf5-p  |zf~p  C*’X(Y,  Z)dXdYdZ. 

Jx> o  ' Y>0  J Z>0  9 

Making  the  transformations  U  =  X  +  Y  +  Z,  R=U  _JsYU ~iit  T  »  if^Zlf**  and 


using  (2.9)  gives 


F  *  Tm(a  +  b+  c,  $) A, 


where  A  is  the  left  hand  side  of  (2.14).  On  the  other  hand,  it  may  be  shown 


r  -  r  (a)r  (b,  *>r  <c.  x)e^xc  fi). 


Hence  (2.14)  follows. 


(ii)  The  integral  is  given  by  the  coefficient  of  C*'*(U,  V)/k!fc!C  (I)  in 


f  /*  a-p 
J  j  R j  etr(ARA’H'UH)dR  etr(aH'VH)dh 

' 0  (m)  J 0 

•  +  *Hl*l*  C.fo  h.x'-’K  lo<m)VAM'U’m> 

C. (aH'Vit)dH/k! I ! (a  +  p)  ,  because  of  Constantine  (1963,  Eq.  (60)), 

A  K 

leading  to  the  desired  result  (2.15)  on  using  (2.2). 


etr(ARA'H'UH)dR  etr(aH'VH)dh 


JK.  - 


3.  JOINT  DISTRIBUTIONS  OF  THE  ROOTS  OF  D  AND  A 

The  latent  roots  of  A  =  (S j  +  Sj)  ^*S  i  < S ^  +  S?)  ^  may  be  utilised 

as  suitable  preliminary  test  statistics  for  the  hypothesis  Ei  "  Z 2» 

where  the  S.  are  independently  distributed  as  Wishart  W  (n.,E.),  i  **  1,  2. 
i  nil 

The  distributions  of  the  roots  of  Sp.  SjS2  have  been  considered  for  the 

K  X 

case  fij  “  0  by  Davis  (1979b)  who  expressed  them  in  terms  of  the  C  *  .  We 

9 

-If  JU 

consider  the  matrix  D  =  (Si  +  S2)  So(Si  +  S2)  in  connection  with 
the  discriminant  problem  under  the  heterogeneity  of  covariance  matrices, 
where  Sq  is  distributed  as  noncentral  Wishart  W^(nQ,  Eq,  fig)  independently 
of  Si  and  S2.  When  Ei  =  E2,  D  and  A  are  independent  and  D  is  the  well- 
known  test  matrix  for  the  null  hypothesis  Ho  =  0  in  MANOVA.  When 
II  +  Z2 »  they  are  dependent  and  we  may  be  interested  in  the  conditional 
distributions  of  the  roots  of  D  or  suitable  functions  of  them  (conditional 
on  the  roots  of  A) , 

Let  di,  d2,...,dm  and  a^  ,  a2,...,am  denote  the  latent  roots  of  D  and 

A  respectively  (d,  >  d  >...>d„  >  0,  1  >  a,  >  a^>-‘->a  >  0).  In  this 

i  /  in  i  z  m 

section,  we  shall  derive  the  joint  density  function  of  di,...,dm  and 
,a  and  the  joint  distribution  functions  of  di  and  aj.  The 
'marginal'  and  conditional  distributions  of  d 1 , . . . ,dm  will  be  considered 

in  the  next  section.  It  is  shown  that  these  distributions  can  be  expressed 

•c  X 

in  terms  of  the  C  ’  for  the  null  case,  Dn  «  0. 

9 

Make  the  transformations 


jl. 

ic  .  c  * 


D  “  S  *S0S  2,  A  -  S  *>  Si  +  S2 


(3.1) 
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in  the  joint  density  function  of  Sq,  Sj  and s 2*  then  the  Jacobian 

,  ,  2P  ^ 

of  the  above  transformation  is  |S|  .  D  and  A  have  the  same  roots  as 

D  ■  h'DH  and  A  =<  H*  AH,  H  e  0(m),  respectively,’  averaging  over  0(m), 
we  have  the  joint  density  function  of  D,  A  and  S  in  the  form,  using 


(2.2), 


-  2  -1  -  ^O-p  -  *2n i“P  ..  Sjnj-p  _1 

f(D,  A,  S)  =  [  JIQ  rm(Ban±)  |2Z±|  ]  |d)  |a|  |i  -  a|  etr(-h  r2  S) 


^4  =  0  ~P  „O0  f  „  \  -  i  -1  -1  If  1  “  f 

|S|  l  (-1)  C*»  (E0  S,  (Ej  -  Z2  )S)C*pA(D,  A)/kUt2rC.(I).  (3.2) 

k,a;<()  <+>  t  9 

Integrating  (3.2)  over  S  >  0,  using  (2.9),  yields  tne  joint  density 
function  of  D  and  A  in  tne  form 


-1  -1  -1 


f(D,  A>  -  C41|D|Ji,,0-p|A|!ini"P|l  -  Ai'in2'11C,X;/“1)f(!iJoni)( 

1  1  «<  •> 

C*’A(£o  £2.  Z?T.2  -  A)/k!t!C(Ji(I), 


there 


C“*  '  rmft1'0V>il!o 

Here  we  snould  acknowledge  tne  fact  that  (3,3)  will  not  Vo  convergent  for 
all  D,  vwhicn  is  a  multivariate  generalization  of  F. 

The  Joint  Density  Function  of  and 


(3.3) 


(3. A) 


,a  • 
’  m 


By  the  usual  method,  the  joint  density  function  of  dj,,..,d^  and 

ai,...,a  is  derived  from  (3.3)  in  the  form 
m 

*1 . «„)  -  c42|D|la,»~p|A|'5ni'p|i  -Al^-^nJ  (dj  -dj) 


»'<*,  -  »,)  C  1  ”.>4e 

£<j  1  J  1=1  i  V  < 


K,A  K,X 


-1 

(L0  l 


2* 


-1 

£l  £2  "  I) 
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(3.5) 


CK(D)Cx(A)/k!£!CK(I)Cx(I), 

i  2  2 

where  C42  “  /r  (*sm)]  with  C41  given  by  (3.4). 

The  Joint  Distribution  Function  of  dj  and  aj . 

We  shall  obtain  the  joint  distribution  of  the  largest  roots  dj 
and  of  D  and  A  respectively,  P(0<dj<6,  0<aj<a)  =  P(0<D<6I,  0<A<aI). 
It  may  be  shown  that 

[5I  |5|Jsno-Pc<.^(D  /^dn  =,  r  (p)r  (!5n0,K)o^xc.(i).,'5mno+k 

Jg  mm  <$>  9  0 


C,(A)/r  (*5no+P«  k)C.(I),  because  of  (2.15),  and  then  (2.4), 
A  m  a 


(3.6) 


lAj'-sni-Pjj  _  A|i?l2"pC,(A)d4  «  P  (p) J  J  T  Ojnj,  t) 

(  0  1  ’  1  1  X  m  ‘’ojteo'X  m 

2  i 

(-ljn-2+p)  (9^*X)  C  d)^1+t/r  (^i+p*  t)s!  ,  because  of  (2.11).  (3.7) 


a  t  t 


Thus,  (3.6)  and  (3.7)  with  (3.3)  establish  the  joint  distribution  func¬ 


tion  of  dj  and  in  the  form 


P(0<d!<«,  0<a,<ct)  -  C.36i,nn”(i'3,ni  IK>J.,(-l)t(lS1E0  14) 


eT'£2  -  n tt»,)t 


-1  -1 


(-4n2+p)„g  (’•  (I)«  /kilieK^no+p)  (>ini+p)_C.  (I) , 

C7(T,XT  r  TA 


where 


(3.8) 


C43  -[  rm<P)l?rm<!S^o  ni>f  rm<!a™2)  "Q  rm(^ni  +  P>lEiZ21|J5niJ_1  • 
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4.  MARGINAL  AND  CONDITIONAL  DISTRIBUTIONS  OF  THE  ROOTS  OF  D. 

Following  the  notation  in  Section  3,  we  shall  derive  the  ’marginal' 

joint  density  function  of  d,  d  ,  the  distribution  function  of  d, 

m 

and  the  density  function  of  t.r  D,  in  this  section.  The  conditional 
joint  density  function  of  ^1,...,^  given  ai,»..,am  and  the  conditional 
distribution  function  P(0<dj<5 | 0<ai<a)  are  obtained  from  these 
’marginal'  distributions  and  the  results  of  Section  3. 

The  joint  density  function  of  D  and  A,  (3.3),  is  integrated  over 
0  <  A  <  I  to  obtain  the  density  function  of  D;  making  the  transformation 
A  -►  H’AH,  II  e  0(m)  ,  and  averaging  over  0(m)  with  the  use  of  (2.10)  and 
(2.4),  and  then  averaging  over  0  <  A  <  I  with  the  use  of  Constantine 
(1963,  Eq.  (22)),  we.  can  establish  the  density  function  of  D  in  the  form 

*  X  -1  -1 

c.  ’(I0  l2,  El  F2  -  DC  (D)/k»«!(,5(n1  +n?)).C  (I),  (4.1) 

Y  K  A  K 

where  C51  =  Tjh  l  n  )  [  r  (Jjn0)r  ftOi!  +  n2))  JJ  |l  Z^l^i]"1. 
i~n  J.^0 

The  Joint  Density  Function  of  d,  d  . 

_  m 

By  the  usual  method,  the  joint  density  function  of  d,  ...  d  is 

m 

readily  derived  from  (4.1)  in  the  form 

f(di,...,dm)  =  C5?|nrl5n°"P  II?  (d.  -  d.)f  (-l)f(Jsni)x(%  l  n  ) 

i<j  J  k,X;4i  1=0  $ 

9^’V^A0:o^2.  eT1E2  -  I)C^(D)/k!  £.!  (Jj(nj  +  njD^U),  (4.2) 


JO 


where 

c52  -  Cs^^/r^Jsm).  (4.3) 

The  Distribution  Function  of  dj. 

The  distribution  function  of  ,  P(dj  <  6)  =  P(D  <  61),  is 
obtained  from  (4.1)  in  the  form,  using  Constantine  (1963,  Eq.  (22)), 

P<dl  <  5>  *  •'i'/*’' 

El ^2  -  I)«k/k!l!fls»o  +  9)Kih<n\  +  n2)) A,  (4.4) 

where 

2  j  j  , 

c 53  -  r  (p)r  Oj  L  n  )[  r  (%n0  +  p)r  (5s(ni  +  n2))  Ii  |E.E2  (4.5) 

m  m  £_o  1  m  1  i^o 

The  Density  Function  of  tr  D. 

By  the  similar  method  to  that  used  for  deriving  the  Hotelling  trace  in 
Davis  (1979b,  Section  8),  the  density  function  of  d  =  tr  D  equal  to 
tr  D,  is  obtained  from  (4.1)  in  the  form 

f(d)  -  C54^n0-r  (“l)f (*sn0)  (*sni)  (%  Z  n.). 

»c ,  A ;  <#>  K  A  i«0  1  * 

e'/c^W^,  Z~\'z2  -  I)dk/kU!(!sinn0)k(-s(nl  +  n2))x,  (4.6) 

2  1  ,  -l.U,  -1 

where  C54  ■  F  ft  I  n  ){  r(*sim»o)r  (‘^(nj  +  n2))  II  |£.£2  I  il 
m  i=0  1  m  i-0 

(4.6)  is  a  generalization  of  Constantine's  series  for  the  Hotelling  trace 
(see  Constantine  (1966)).  Presumably  it  may  converge  only  for  |d | <1 . 

The  Conditional  Distributions  of  the  Roots  of  P  Given  Those  of  A 

It  is  easily  shown  that  the  joint  density  function  of  the  roots  S),...,a 
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of  A  is  given  by 


f(alt...,a  )  =  C55|A|isni‘p|l  -  Al^-P  nT(a  -  a  ) 
m  i<j  J 

-1 

lF0(!s(ni+n2);  Z2  -  I,  ~  A), 

1  2  2  1  1 

where  C55  -  r  (*s(m  +  nz))^  [  r  (*jm)  n  r  (ijn  > - 1 r  1E2  i^1]"  ,  and 
m  in  i~l 

that  the  distribution  function  of  aj  is  given  by 


(4-7) 


P(aj  <  a)  ■  C56  [o(l  -  a)  ]'£mn l 2f  1  Ojn  1 ,  Js(n1+n2);  %ni+p;  -a(l-o)  I^1!^).  f  .8) 

where  Cs6  -  rm(p) ^(-^(ni-H^)  )  Irm(!Jn2) ^(^l  +  P)|  £1^2  I  5”1]  ,  and  the 

2F1  is  the  Gaussian  hypergeometric  function  of  matrix  argument  (see  e.g., 

James  (1964)). 

Hence,  the  conditional  joint  density  function  of  d-j,.,.,  d^  given 

ai,...,a_  is  obtained  as 
*  m 

f(di,...,dm|a1,...,am)  -  f(d1,...,dm,  al»* • • »aJ /f <al •  • • • *ara> .  (A,9) 

where  f(dlt...,dm,  aj,...^)  and  f(»i . a^)  are  given  by  (3.5)  and 

(4.7)  respectively.  When  Ej  =  E2,  (4.9)  becomes  independent  of  A  in  the 
form  of  the  density  function  of  the  roots  of  the  F  matrix,  expressed  in 
terms  of  the  jFg  hypergeometric  function  with  two  matrix  arguments  (see 
e.g.,  James  (1964.,  Eq.  (65))). 

The  conditional  distribution  function  of  dj  given  aj  is  obtained  as 
P(dj  <  d|«j  <  a)  -  P(dj  <  5,  aj  <  o)/P(aI  <  a),  (4.10) 

where  P(dj  <  $,  Sj  <  a)  and  P(ai  <  a)  are  given  by  (3.8)  and  (4,8) 
respectively.  When  Ei  *=  E2,  (4.10)  reduces  to  the  distribution  function 
of  the  largest  root  of  the  F  matrix,  expressed  in  terms  of  the  jFi  hyper- 
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geometric  function  of  matrix  arguments  (see  e.g.,  Chikuse  (1977, 

Eq.  (2.4))). 

Before  closing  this  paper  it  is  noted  that  we  can  extend  some 
of  these  results  by  introducing  invariant  polynomials  with  larger 
numbers  of  argument  matrices,  extending  the  work  of  Davis  (1979a), 
(1979b).  These  will  be  discussed  in  subsequent  papers  (see  e.g., 
Chikuse  (1980)). 
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